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Abstract. We study theoretically the class of mixed-valent Kondo insulators, 
employing a recently developed local moment approach to heavy Fermion systems 
using the asymmetric periodic Anderson model (PAM). Novel features in spectra 
and transport, observable experimentally but lying outside the scope of the 

Q^ ' symmetric PAM or the Kondo lattice model, emerge naturally within the present 

theory. We argue in particular that a shoulder-like feature in the optical 
conductivity, that is distinct from the usual mid-infrared or direct gap peak and 

C^ , has been observed experimentally in mixed-valent compounds such as CeOs4Sbi2 

and YbAl3, is of intrinsic origin. Detailed comparison is made between the 

C^ ' resultant theory and transport/optical experiments on the fiUed-skutterudite 

compound CeOs4Sbi2, and good agreement is obtained. 



PACS numbers: 71.27.-|-a Strongly correlated electron systems; heavy Fermions 
- 75.20.Hr Local moment in compounds and alloys; Kondo effect, valence 
fluctuations, heavy Fermions 



1. Introduction. 



p^ ' Kondo insulator compounds constitute a sub-class of lanthanide/actinide based heavy 

Fermion materials. Characterised by small spectral and optical gaps, and an 

^^ ' activated resistivity at low temperatures, they have been under intense theoretical and 

experimental investigation for many years (for reviews see [1-7]). From a theoretical 
perspective, the essential framework for understanding heavy Fermion materials is the 

t . periodic Anderson model (PAM) , wherein a single correlated /-level in each unit cell 

of the lattice hybridizes locally with a non-interacting conduction band. 

Theories for Kondo insulators (KI) have generally been based on the particle- 
hole symmetric limit of the PAM [1,7-13], in which the /-level occupancy (uf) and 
the conduction band filling (ric) are each equal to unity. Although Uf = 1 = Uc 
C^ ' indeed satisfies the condition for an insulating ground state, it is not the generic case, 

which by contrast is n/ + ric = 2 (as shown and discussed in [3, 10, 14]). The strong 
coupling limit of the PAM is of course the Kondo lattice model with wholly localised 
/-electrons and hence Uf — 1 necessarily, so one naturally expects that if a particular 
KI is sufficiently strongly correlated then the particle-hole symmetric PAM should 
provide a sound description of it. And indeed theories baseed on the symmetric PAM 
have been able to describe, even quantitatively, many aspects of a number of these 
materials (see e.g. [13]). 

One cannot however expect such theories to be complete, since the fact that 
the generic condition for Kondo insulating behaviour is Uf = 2 — ric [3, 10, 14] itself 
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suggests that KI materials are likely as a rule to be mixed- valent [nj ^ 1); as indeed 
seems to be the case experimentally [3,7], and which behaviour lies beyond the scope 
of the symmetric PAM. A more general treatment of KIs to encompass mixed- valent 
behaviour is thus clearly desirable, based on the general asymmetric PAM with the 
constraint n/ = 2 — nc 7^ 1. It is this we consider here. 

In a previous paper [15], working within the general framework of dynamical 
mean-field theory (DMFT) [8,9,16,17], we developed a local moment approach (LMA) 
(see [12, 13, 18-25] for details) to heavy Fermion systems. This encompassed both 
strong coupling Kondo lattice behaviour (n/ ^- 1, ric arbitrary), as well as the mixed- 
valence regime (n/ and Uc arbitrary). The theory was subsequently compared to 
experiments on several heavy Fermion metals, and excellent quantitative agreement 
was found [26]. Although comparison was made only to metallic systems, the 
insulating ground state is just a particular solution of the same basic theory, that 
may be obtained by restricting the parameter set to the line Uf + Uc ~ 2. In the 
present paper we treat the generic case by supplementing the theory of [15] with a 
constraint on the total filling, rif -\- Uc ~ 2, but allowing Uf and Uc to deviate from 
unity. This enables us to capture the mixed-valency along with the insulating nature 
of these systems. 

We begin with a brief overview of the model and formalism used. Theoretical 
results for spectra and transport are discussed in section 3, followed by a detailed 
comparison with experiments on the filled skutterudite compound CeOsiShi2, and a 
brief conclusion in section 4. 

2. Model and theory. 

In standard notation, the Hamiltonian for the PAM is given by: 

+ Vj2{flc,,+h.c.) (2.1) 

The first pair of terms describe the non-interacting conduction (c) band. The first 
gives the kinetic energy (or 'free' conduction band H^), with the nearest neighbour 
hopping tij = —t scaled as i^ oc t* /y/Z^ in the large dimensional limit of coordination 
number Zc —* 00 [8,9,16,17], and t* = 1 taken throughout as the unit of energy; while 
the second gives the c-orbital energy, such that varying ec simply shifts the centre of 
gravity of the free conduction band relative to the Fermi level, and as such controls the 
conduction electron filling Uc- The third term in H represents the /-orbital energy (e/) 
and the on-site Coulomb repulsion (U) for the localised /-orbitals, and the final term 
denotes the local hybridization between the c and /-electrons which is responsible for 
making the otherwise localized /-electrons itinerant. 

The free conduction band may be diagonalised, H^! = X^k o- ^k'^ko-'-ko-' with cor- 
responding density of states po(e) — A^~^ X)k^(^ ~ ^k)- While the basic formalism is 
valid for any underlying lattice and associated po{e), we consider in this paper the spe- 
cific case of the hypercubic lattice (which is Bloch-decomposable, unlike e.g. a Bethe 
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lattice), for which po{e) = cxp(— e^)/Y^ is a Gaussian. 

Within DMFT, the /-electron self-energy S/(w;r), representing many-body 
scattering due to electron interactions, is rendered purely local (i.e. site-diagonal or 
momentum independent). The local, site-diagonal Green functions for the conduction- 
(G^(w)) and /-electrons {G^ {oj)) are then given by [8,9, 14, 16, 17] 

G^M = / de— P-^^, EE / depo(e)G=(e;c^) (2.2a) 



G-^(a 



1 






(2.26) 



UJ+ -Ef -j:f{u;;T) 

(with Lu^ — Lu + iO^). Solution of these DMFT equations naturally requires 
a knowledge - and self-consistent determination - of the /-electron self-energy 
T,f{uj;T), which DMFT by itself does not of course prescribe. For this, we employ 
the physically transparent local moment approach (LMA) [12-15, 18-26]. This 
handles non-perturbatively all interaction strengths from weak to strong coupling, 
and all relevant energy/temperature scales, while at the same time recovering the 
dictates of Fermi liquid behaviour at low energy/temperature scales. Originally 
introduced to describe Anderson impurity models [18-25], for which it is found to 
agree well with e.g. numerical renormalization group calculations and a number of 
exact results, lattice-based heavy Fermion systems have also been considered within 
DMFT-FLMA [12-15,26]. 

Here we sketch in brief only the essential elements of the LMA: 
(i)The starting point is simple static mean-field, i.e. unrestricted Hartree-Fock. This 
has the virtue of recognising local moment formation as the first effect of electron 
interactions, by introducing the possibility of local moments from the outset. But by 
itself it is inadequate, for two reasons. First, it results in a (locally doubly-degenerate) 
symmetry broken mean-field state, which is not perturbatively connected (in U) to 
the non-interacting limit and in consequence violates Fermi liquid behaviour at low 
energies. Second, its inherently static nature cannot by construction capture electron 
correlation effects. It is these signal limitations the LMA overcomes, 
(ii) Electron correlations, embodied in dynamical self-energies, are incorporated within 
the framework of a spin-rotationally invariant two-self-energy description which is 
an inevitable consequence of the underlying two mean-field saddle points (and from 
which the conventional single self-energy S/(w; T) is recovered simply as a byproduct). 
The resultant dynamical self-energies are built diagrammatically from, and are self- 
consistently determined functionals of, the underlying mean- field propagators. They 
include in particular a non-perturbative class of diagrams that capture the spin-flip 
dynamics central to the physics of the PAM. 

(iii)The third, key element of the LMA is that of symmetry restoration [12, 19,20,24]: 
self-consistent restoration of the symmetry broken at pure mean-field level, and hence 
correct recovery of the low-energy local Fermi liquid behaviour that reflects adiabatic 
continuity to the non- interacting limit. This is embodied in a single self-consistency 
condition on the two-self-energies precisely at the Fermi level (w = 0), which in practice 
amounts to a self-consistent determination of the local moment (supplanting the simple 
'gap equation' for such that arises at crude mean-field level). 

(iv) Luttinger's theorem [2,27], itself a reflection of perturbative continuity to the 
non-interacting limit (i.e. Il — liaij_^duj {dY^f{uj;T — Q)/duj)G^ {cj) — 0), is 
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also satisfied by the LMA [14]. We note too that adiabatic continuity to the non- 
interacting limit is characteristic of both the heavy Fermion (metallic) and Kondo 
insulating states [12,13]. 

Full details regarding the structure and implementation of the Local Moment 
Approach for the generic asymmetric PAM are given in [14, 15], to which the reader 
is referred for further information. 

Since our objective is to study generic Kondo insulators, corresponding to 
ric + n/ = 2 (for T = where the distinction between an insulator and a metal has 
strict meaning), we now examine the general conditions under which an insulating 
ground state is obtained. To this end we note that the Luttinger theorem II — can 
be expressed as the following exact statement for the Fermi surface of the PAM, as 
shown in [14]: 

i(n, + n/)=/" " ' po{e)de + 9{-e}) (2.3) 

■.Z — oo 

(with 9{x) the unit step function). Here e*f — eZ/V'^, and eZ- is the renormalised or 
effective /-level given by e) = ef + T,f{uj = 0; T = 0) (where S^ = ReS/). 

From this it follows that, for any conduction band poi^), the KI filling constraint 
Uc + Uf = 2 is satisfied ii e*f = (whether eZ = 0+ or 0^). Moreover for a non-compact 
bare density of states po{e), such as the Gaussian appropriate to the hypercubic lattice 
we consider, e j- = is the only possibility that satisfies Uc + Uf = 2. As such, e*f = 
is the general condition for a KI that we seek; with the (particle-hole) symmetric KI 
(e/ = —U/2 and ec = 0) simply the particular case for which ric = 1 = n/. 

The condition el — also leads naturally to a gap at the Fermi level {ui = 0) in 
the T = single-particle spectra D'^{ui) and D^{uj). This is most easily seen from the 
limiting low-frequency behaviour of equations (2.2a) and (2.26) for T ~ 0, obtained 
from a \ow-lu quasiparticle expansion [14] and given by 

D^Lu) ^ poi-e, - ^^) (2.4a) 

V'Dficj) ^ \^., Po{-ec ^^) (2.46) 

[ll) — e^y u! — ef 

where uj = lo/ujl. Here, ujl = ZV'^ /t^ is the characteristic low-energy Fermi liquid 
scale in the problem (with Z = [1 — (9E?(a;; T ~ 0)/9w)^=o]~^ the usual quasiparticle 
weight, or inverse mass renormalization factor). From this it is readily seen that with 
e*f = 0, D'^{uj = 0) = (and likewise for the /-spectrum), as one expects for an 
insulator. 

In this paper we implement the filling condition ric + Uf = 2 simply by 
supplementing the LMA with the T = constraint e* = e/ -I- i;^(w = 0;T = 0) = 
(which is algorithmically simple), ensuring thereby direct access to the generic Kondo 
insulating states of interest. The LMA yields directly the local Green functions (G"^ 
and G^) and self-energies, knowledge of which is well known to be sufficient within 
DMFT [8,9,16,17] to determine d.c. transport and optical properties, as detailed e.g. 
in [13,15]. In the next section, we discuss our theoretical results. 
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Figure 1. Variation of /-level asymmetry rj with |ec| (= — £c) for the Kondo 
insulating state, for fixed V^ = 0.2 and three different interactions U = 2.6 (solid 
line), 3.8 (dashed line) and 5.2 (point-dash). With increasing jed, the insulator 
becomes progressively more mixed-valent (the /-electron filling dropping from 
rif = 1 for \tc\ = to rif ~ 0.63 for \ec\ = 0.4). 



3. Theoretical results and discussion 



The PAM, equation 2.1, is characterised by four 'bare' material parameters, U,V,ec 
and the /-level asymmetry r] = 1 + 2e//C/ (or equivalently ej itself). The local spectra 
naturally depend on these parameters, and the (T = 0) total filling is given by 



Uf 



duj [D'{uj; U, V, e,, r,) + D^ {oj- U, V, e„ r])] 



(3.1) 



For a metallic phase, (U,V,ec,r]) are in general independent parameters. But for the 
Kondo insulating state, the constraint ric + Uf = 2 obviously implies that only 3 of 
the bare parameters are independent, e.g. 77 = r](lJ, V, Ec)- 

We illustrate this is figure 1, showing the resultant rj vs jed (with Cc < and 
hence Uc > 1), for a fixed value of the hybridization V and three different interactions 
U. All curves meet at the origin, this being the symmetric KI (n/ = 1 = Uc) for 
which 77 = = ec for all U and F, i.e. 77(^7, V, ec = 0) = is independent of U and 
V . Away from Cc — Q however, 77 depends generically on all of U, V and Ec. On 
increasing jed for any given U and F, as in figure 1, the system becomes progressively 
mixed-valent: 77 increases (i.e. the /-level moves upwards towards the Fermi level), 
and rif steadily diminishes from tt./ = 1 at ec = 0, with a concomitant increase in 
the resultant quasiparticle weight/inverse mass renormalization Z (and hence in the 
low-energy scale lol = ZV"^). In fact over the [/-range shown in figure 1, the resultant 
n/s for given Cc barely change with C/, Uf dropping to ~ 0.63 for led = 0.4. 

To illustrate typical local single-particle dynamics, figure 2 shows the T = LMA 
c-and /-spectra vs Cj = uj /ujl (left and right panels respectively) for a fixed Cc = —0.4, 
U = 2 and a range of different hybridization couplings V as indicated. The associated 
Hf differs only slightly with V over the range shown, with Uf ~ 0.6 indicating mixed- 
valent character. The single-particle spectra for this representative asymmetric KI are 
indeed seen to be gapped at the Fermi level, with magnitude ~ wl, and the obvious 
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key point is that the gap is asymmetrie about the Fermi level (it is of course strictly 
soft for a hypercubic lattice, but with exponentially small spectral density in the 
vicinity of the Fermi level). The four sets of spectra are seen to be quite distinct away 
from the Fermi level, as one expects. However the limiting low-frequency spectral 



forms are given (asymptotically exactly) by equations (2.4a, 2.46) which, with e 



/ 







as appropriate to the KI, are seen to depend solely on tc and Co; i.e. they should be 
independent oiVoiU (which we note has nothing per se to do with the 'universal 
scaling' of spectra characteristic of the strong coupling limit n/ ~ 1 [12-15]). That 
the LMA correctly recovers this behaviour correctly is self-evident from the figure. 
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Figure 2. The local conduction electron spectra (left panel) and /-electron 
spectra (right panel, with the constant Aq = nV^po{—ec)) are shown vs. 
Gj = uj/u)L for a fixed U = 2.0, ec = -0.4 and varying V^ = 0.005 (solid line), 0.01 
(dotted), 0.03 (dashed) 0.075 (double dot dashed). An asymmetric gap is seen 
straddling the Fermi level. Inset: magnified view of the low frequency /-spectrum. 



3.1. Optical conductivity. 

We consider now the frequency dependence of the optical conductivity a{uj]T)^ 
focusing specifically on T = (results at finite-T will be included in the following 
section). Figure 3 shows the optical conductivity obtained from the LMA for a 
representative mixed- valent Kondo insulator. The dominant feature in the optics is the 
usual strong, direct gap absorption ('mid-infrared peak'). In the absence of scattering 
due to electron interactions there would be no absorption whatever below the direct 
gap [13, 15]; but, just as for the (particle-hole) symmetric KI considered in [13], 
interaction-induced many-body scattering leads as seen to continuous absorption all 
the way down to the indirect gap scale to = Ajnd = Sojl (which is why 17(0'; 0) is shown 
as a function of a) = '^I'-Ol-, figure 3 showing there is indeed negligible absorption below 
a;~2). 
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Figure 3. Le/4 'panel: Optical conductivity a-(a;; 0) ws tj = Lj/ajj^ for a 
representative rnixed-valent Kondo insulator (ny ~ 0.73), with IJ = 3.8, V^ = 
0.2, ec = —0.3 and r] = 0.8. Right panel: the renormalized bandstructure <.i;(ek) 
lis Ek (solid lines). For the dotted lines, see discussion in text. 



Figure 3 also shows that the optical conductivity possesses a distinct low- 
frequency shoulder, lying somewhat above the indirect gap (at J; « 10 in the present 
example). Such a feature occurs neither in the non-interacting limit nor for the 
interacting symmetric Kl [13]. From investigation of a wide range of underlying 
material parameters however, we find it to be entirely typical of interacting mixed- 
valent (necessarily asymmetric) Kondo insulators, albeit naturally becoming less 
pronounced as the symmetric limit is approached. In previous work [15, 26] we 
have also found such behaviour to be characteristic of correlated intermediate valence 
metals. We conclude therefore that a low-frequency shoulder should typically exist, as 
an intrinsic optical feature, in interacting mixed-valent systems whether metallic or 
insulating; and that it arises from a combination of many-body scattering and mixed- 
valency. Experimentally, such a shoulder has been observed in the optical conductivity 
of Kondo insulators such as CeOs4Sbi2 [28] as well as in intermediate valence metals 
such as YbAla [29] . In [26] we considered the case of YbAla in detail, and obtained 
very good agreement between theory and experiment, including striking reproduction 
of the low-energy optical shoulder. Analogous comparison for CeOs4Sbi2 will be given 
in the following section. 

First, however, we seek to interrogate the optics in more microscopic detail. In 
the absence of vertex corrections (as appropriate to DMFT [8,9, 16, 17]) the T = 
optical conductivity is given by [13, 15] 

(3.2a) 



(3.2&) 



ct(w;0) cx / de po{e) /(e;cj) 

(bar extraneous constants), with 

1 /•" 
I{e;uj) ^ - / duji D%e; 101)0" {e;uji+uj). 

Here, D'^(e;uj) = D'^(ek — ^',^) is the ek-resolved conduction electron spectrum 
(accessible e.g. via ARPES): D'^{e;uj) — — (l/7r)ImG"^(e;u;), with the k-space 
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propagator C^ {e\^-, uj) given by 

G\e^-u)^[-i{oj)-e^]-^ (3.3a) 

where 

^(uj) = LO'^ — £„ \ -^T—, 7^ (3.3&) 

(such that the local, site-diagonal propagator G'^(cj) = A^^^ ^j, G"^(ek; tj) = 
j de po{e)G'^ {e; uj) as in equation (2.2a)). 

The physical content of equations (3.2) is clear: for given e = e^ (as in 
equation( 3.26)), optical transitions occur from a state below the Fermi level (lying at 
Wi < 0) to a state oj higher in energy that hes above the Fermi level (at oji+uj > 0); all 
transitions being 'direct' in the sense that absorption occurs for given/fixed e^. The 
net optical conductivity at frequency lu is then, as in equation (3.2a), the sum of all 
such processes over the full range of e = ek (with density po(e) — N^^ J2k <^(^k — e))- 
This picture is of course wholly familiar at an elementary level in the context of 
non-interacting electrons (or 'renormalized' non-interacting electrons, as considered 
below). But we emphasise that it is quite general: it applies equally well to the 
fully interacting case, where the states between which absorptive transitions occur are 
many-body states. There is of course an obvious difference between the interacting and 
non-interacting cases, namely the existence of scattering (and hence 'lifetime' effects) 
due to electron interactions in the former case. That in turn generates a profound 
difference between optics in the two cases, which we now consider since it also throws 
light on the optical conductivity shown in figure 3. 

The simplest description of the optics is at the level of 'renormalized 
bandstructure' [9]. Here, 7(0;) in equation (3.3a) is taken to be purely real, 7(0;) = 
7^(w), neglect of 7^(0;) — Im7(u;) meaning from equation (3.3&) that the imaginary 
part of the /-electron self-energy — the source of all scattering — is completely 
neglected. In addition, 7^(w) is further approximated by its asymptotic low- frequency 
behaviour, which (from equation (3.36)) comes from that for S?(cj;0) = ReS/(w;0), 
namely the simple Taylor expansion S?(a;;0) ^ S?(0;0) — (1/Z — 1) w (with Z the 
usual quasiparticle weight). With this approximation G"^(ek;w) ~ [7^(0;) — ek]~^ has 
poles at two frequencies, lo — cij+(ek) > above the Fermi level and w = io~(e.\^ < 
below it. These are given explicitly by 



U! 



(ek) = | (ee + ek) ± Vi^c + ek)^ + 4ZV^ (3.4a) 



(where we use the fact that the renormalized level e*f =^ ef + S^(0; 0) = for a generic 

KI, as discussed in section 2), with a gap between them A(ek) = cij+(ek) — "^ (^k) of 

A(ek) = V(ec + ek)2+4Zy2. (3.45) 

Equation (3.4a)) is simply the renormalized bandstructure of the problem (reducing 
trivially to the non-interacting limit result for Z = 1), the two branches Lo~^{eu) and 
'^"(sk) reflecting physically the fact that, per unit cell, a single /-level hybridizes 
locally to a single conduction band. The minimum gap between them, the 'direct gap' 
at this level, occurs for ek — —t^ and is thus A^jr = A(ek = — e^) ~ 2\fZY . 

Since the resultant Z)'^(ek;tj') ex ImG'^(ek;ti;') contains two poles, at w' — w+(ek) 
and a;^(ek), it follows from equation (3.26) that optical absorption can only arise from 
transitions between these two (5-peaks — provided of course their separation (^"""(ek) — 
|^~(ek) concides with the requisite absorption frequency w; and in consequence that 
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Figure 4. For the same bare parameters as in figure 3. Left 'panel: T = 
conduction electron spectrum D'^(e]^\LD) vs ui, for fixed e^ = 0.3 (= — ec)- The 
solid arrows indicate the location of the pure poles that would arise from a simple 
rcnornialized band picture (following from the renormalizcd bandstructure shown 
in figure 3, right panel). An additional spectral feature above the Fermi level, 
arising from interactions, is indicated by a single arrow. Right panel: The same, 
shown for an e^ = 0.6. 



there is no absorption at all (i.e. a{(jj; 0) — 0) for uj < Adir — 2\fZV . This level of 
description may be improved, very slightly, by retaining C^ies^^uS) ~ Yi^{^) ^ ^k]^^ 
but with the full LMA 7^(w) employed instead of its asymptotic low-w expansion. 
The results of such a calculation are shown in figure 3 (right panel, solid lines, where 
the 'band branching' seen over a narrow range around ek = is simply a consequence 
of retaining 7^(w) but neglecting 7^(0;)). 

The above 'renormalizcd band' description, commonly employed though it is, 
is qualitatively inadequate. The picture changes markedly when interaction-induced 
scattering is properly retained, as embodied in the imaginary part of the /-electron 
self-energy and hence a non- vanishing conduction electron scattering rate 7^(w). To 
illustrate this, figure 4 shows the full LMA conduction electron spectra D'^{e\^\i^) vs 
Lu for two values of e = Ck, viz ek = — Ec = 0.3 (left panel) and ek = 0.6 (right panel). 
At the simplistic level of renormalized bands, as in the right panel to figure 3, each of 
these spectra would consist of one pole below the Fermi level (w = 0) and one above it. 
Reality is clearly different: the spectra are significantly broadened due to interactions 
and form a continuum (save for the expected gap in the immediate vicinity of the 
Fermi level). It is for this reason that, in contrast to the renormalized band picture, 
conduction electron spectra for essentially any ek contribute to optical absorption for 
all frequencies down to the lowest energy (indirect) gap scale. 

In figure 4 we also mark (by solid arrows) the locations in Z)'^(ek;tj) of the two 
nominal poles that would arise at the renormalized band level (their positions can 
be read off from the right panel in figure 3 at the appropriate ek). These are seen 
to correspond rather accurately to the position of the dominant peak maxima in 
Z)''(ek; to). Just above the Fermi level, an additional small peak (marked by an arrow) 
can also be seen in figure 4 . This feature does not of course have any counterpart at 
the renormalized band level, and we find its existence to be characteristic of mixed- 
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Figure 5. For the same bare parameters as in figure 3. Left panel: the T = 
optical conductivity (T{ijj; 0) vs lu (solid line) and the contribution to it arising 
from e = e^ > 0.1 (point-dash), < ek < 0-1 (dashed) and ek < (dotted). Right 
panel: showing further the contributions to (7(uj; 0) from 0.1 < tk < 0-5 (dashed 
line) and ek > 0.5 (dotted). Inset: shown on a lowcr-a; scale. 



valent Kondo insulators (it does not occur in the particle-hole symmetric limit). These 
spectral features, two peaks above and one below the Fermi level, are found to be 
characteristic of the conduction spectra D'^{ei^; ui) for all e^. The positions of the two 
peak maxima above the Fermi level in the full I?'^(ek;w) can be mapped out as a 
function of ek, and are shown in figure 3 (right panel, dotted lines). Note in particular 
that the lower-energy such peak lies close to, but slightly above, the Fermi level for 
all Ck- We thus expect (see equation (3.2&)) significant optical absorption to this peak 
from the dominant spectral peak in D'^{ek;Lo) below the Fermi level; particularly at 
relatively low frequencies where (see figure 3 right panel) for Ck ^ 0.3 or so the latter 
peak itself lies close to the Fermi level, exhibits only modest dispersion with ek, and 
is spectrally sharp (as in figure 4 right panel). 

We return now to the optical conductivity a{uj; 0) given by equations (3.2a, 3.2&), 
shown in the left panel of figure 3 for the representative bare parameters indicated, and 
again on a linear-a> scale in figure 5 (solid line). With the preceding comments in mind, 
our aim is to determine what range of values of e = ek give the primary contribution 
to (j{lo; 0) in different frequency intervals — in particular, frequencies in the vicinity 
of (a) the dominant direct gap absorption (lu w 0.4 in the present example), (b) the 
low-frequency shoulder (lu « 0.1), and (c) the lowest frequency scales down to the 
indirect gap {lu « 0.02). 

To that end, we simply partition the e-integral in equation (3.2a) for cr(Lu; 0), into 
contributions from different e-intervals. This is shown in the left panel to figure 5, 
where the separate contributions to the total ct(w; 0) (solid line) from e = ek > 0.1 
(point-dash line), < ek < 0.1 (dashed) and ek < (dotted) are shown. From this 
it is clear that the dominant contribution to ct^lu; 0) in all three oj-regimcs of interest 
arises from ek > 0.1 (the associated renormalized bandstructure being shown in the 
right panel of figure 3). To distinguish the w-regimes, the right panel in figure 5 
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shows the contribution to a{uj]Q) from 0.1 < Ck < 0.5 (dashed Unc) and Ck > 0.5 
(dotted). From the renormaHzed band description discussed above, the direct gap 
occurs at ek = — Cc (= 0.3 in the present example); and consistent both with this 
and the inevitable broadening induced by scattering, the dominant contribution to 
17(0;; 0) in the vicinity of the strong direct gap absorption at lo k, 0.4 is indeed seen to 
arise from the interval 0.1 < Ck < 0.5 — i.e. ek = — Cc ± 0.2. In contrast, as seen in 
particular from the right inset to figure 5, absorption on the lowest frequencies down 
to the order of the indirect gap is controlled by ek ^ 0.5. This in turn is consistent 
with the form of D'^(ek;w) for Ck's in this range, exemplified by the right panel in 
figure 4 and discussed above. Finally, the low-w shoulder in a{w, 0) (at w « 0.1) that 
is typical of the mixed- valent Kondo insulator, is seen to stem mainly from the interval 
0.1 < ek < 0.5, its existence reflecting the small spectral feature in D'^(ek;w) above 
the Fermi level discussed above, and shown in figure 4. 

4. Comparison to transport in CeOs4Sbi2 

The filled-skutterudite compounds RT4X12 (R=rare earth, T=Transition metal and 
X=pnictide) have attracted much experimental interest, due in part to their possible 
applications as advanced thermoelectric materials [30]. Recently, several groups 
have reported investigations of a cerium based filled-skutterudite Kondo insulator, 
CeOs4Sbi2. Here we make comparison of the theory outlined in this paper to 
experimental measurements of d.c. and optical transport in CeOs4Sbi2. 

We consider first the d.c. transport. In this compound the phonon contribution to 
the d.c. resistivity is quite significant, so for comparison to theory we need to extract 
the magnetic contribution, Pmag{T), from the d.c. resistivity measured directly. This 
is achieved by subtracting the phonon contribution - itself estimated conventionally 
as the resistivity of LaOs^Sbu [31, 32] (see [26] for a detailed discussion) - from 
the experimentally measured d.c. resistivity [33]; and is shown as circles in figure 6. 
Another experimental group has also reported Pmag{T) itself [31]; and although their 
absolute magnitudes are quite different from those of [33], a simple y-axis rescaling 
is sufficient to collapse the two sets of data, as indicated by squares and circles in 
figure 6 (and indicating simply distinct sample geometries in the two cases). 

The low temperature (T < 30K) resistivity can be fitted to a exp[{T*/Ty/^] 
form, suggesting the dominant mechanism of transport in this temperature range to 
be Efros-Shklovskii type variable range hopping [34] (which extrinsic behaviour is not 
of course included in the present theory). A shallow maximum is seen in figure 6 at 
T ~ lOOK, beyond which Pmag{T) decays monotonically. A rough estimate of the low 
energy scale (wl = Z^^/t,) in this system may be obtained through an estimation 
of the (indirect) gap in the experimental optical conductivity [28], shown in the left 
panel of figure 7. This yields Ai„d ~ 15meV ~ 175K, which thus gives the low energy 
scale ujl = Aind/2 ~ 90K (see also the left panel of figure 3 and [13]). 

An idea of the parameter regime this system belongs to can be gleaned 
from two qualitative features evident in the (a;,T)-dependence of the experimental 
o'cxpI'-^'tT) [28] shown in figure 7. First, the highest temperature at which (Toxp('-^;T') 
is measured is 295^ which is ^ Slol, while the optical conductivity at the same 
temperature is affected up to ^ 0.3eV ~ 40^;^. The fact that the thermally induced 
spectral weight redistribution is over a range of frequencies much higher than the 
temperature at which aexp(w;T) is measured is characteristic of relatively weak 
correlations [13]. Second, a closer look at the low-frequency structure of a{co; T = 8K) 
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Figure 6. Experimental d.c. resistivity of CeOs4Sbi2 with the phonon 
contribution subtracted (circles [33] and squares [31]) compared to the 
theoretically determined p{T) (solid line), which has been rescaled by u>l = 86K 
and l/cTQ = 74:fiQcra. The inset shows the corresponding d.c. conductivities. 



reveals an additional absorption feature at a; ^ SOmeV, which in the context of the 
theoretical results obtained above (see e.g. figure 3), is suggestive of mixed- valence 
character. 

With this in mind, wc choose an ec — —0.3, and for various U, V^ calculate 
a{Lu;T). We find that the d.c. resistivity {ui = 0) and the optical conductivity 
determined at C/ ~ 5.5, V^ = 0.5 match best with the experiment, with ujl ~ 86^. 
With these parameter values the /-level asymmetry rj is found to be 0.9, implying 
Uf = 0.75 and Uc — 1.25 consistent with mixed valent behaviour, while the 
quasiparticle weight Z ~ 0.083 implying a relatively low effective mass ~ 12 consistent 
with modest correlations. The resultant theoretical d.c. resistivity is shown in figure 6 
(solid line, with the a:-axis scaled by cj^ — 86K and the y-axis by I/cto — 74/j,r2c77i) . 
And comparison to experiment is seen to be rather good for T > 30K or so (recall as 
above that transport for T < 30K is dominated by variable range hopping [34] , which 
is naturally absent from the theory). 

The theoretical optical conductivity evaluated for the same model parameters as 
in figure 6 is shown in the right panel of figure 7, for similar temperatures as the 
experiment [28] (left panel). An arrow marks the additional absorption seen as a 
shoulder, which arises at w ~ AOmeV and thus compares well with the experimental 
value of ~ 30meV. The theoretical direct gap peak appears at w ~ lOOmeV, while 
the corresponding experimental peak position is ~ TOmeV. The overall theoretical 
lineshape and its thermal evolution also matches rather well with experiment. We 
add further that this comparison does not depend crucially on Ej, being equal to 
—0.3. Varying Ej, by ±0.1 does not change the qualitative picture, although optimal 
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Figure 7. The left panel shows the experimentally measured optical conductivity 
of CeOs4Sbi2 [28] for various temperatures. The right panel shows the 
theoretically determined cf{ll>;T) for the same model parameters as in figure 6 
and similar temperatures to experiment. 



quantitative agreement arises for Cc = —0.3, U = 5.5 and V"^ = 0.5 as employed above. 



5. Conclusion 

We have described in this paper a many-body theory for mixed-valent Kondo 
insulators, employing a local moment approach to the periodic Anderson model within 
the framework of dynamical mean field theory. Kondo insulators were argued as a 
rule to be mixed-valent, with n/ = 2 — ric 7^ 1, and which regime of behaviour is not 
captured by the particle-hole symmetric limit of the PAM (n/ = 1 = Uc)- To that 
end we have considered the general asymmetric PAM, together with the constraint 
Uf -\- Uc = 2 which ensures an insulating gap in the (T = 0) single-particle spectrum 
and related dynamics. 

Single-particle spectra, as well as optical and d.c. conductivities, have been con- 
sidered, and exhibit features specific to mixed-valent behaviour. The k-resolved con- 
duction electron spectra are found to contain additional absorption features above the 
Fermi level, giving rise in turn to an intrinsic shoulder-like absorption in the optical 
conductivity at frequencies lower than the direct (mid-infrared) gap. We emphasise 
that such a structure is found to be characteristic of the interacting mixed-valent 
case, and not restricted to either an insulating or metallic ground state. Thus we be- 
lieve that the shoulder-like feature seen experimentally in the optical conductivity of 
mixed-valent materials such as the insulator CeOs4Sbi2 and the metal YbAla, arises 
intrinsically from a combination of many-body scattering and intermediate valence. A 
direct comparison between theory and transport/optical experiments on CcOs4Sbi2 
has been given. This yields good quantitative agreement, both reaffirming the view 
of CeOs4Sbi2 as a mixed-valent hybridization gap material [33,34] and showing that 
the theory described here for generic mixed-valent Kondo insulators can account for 
the transport/optical properties of these systems. 
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